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Abstract. The flow behavior of a semiflexible polymer in microchannels 
is studied using Multipartiele Collision Dynamics (MFC), a particle- based 
hydrodynamic simulation technique. Conformations, distributions, and radial 
cross-streamline migration are investigated for various bending rigidities, with 
persistence lengths Lp in the range 0.5 < Lp/Lr < 30. The flow behavior is 
governed by the competition between a hydrodynamic lift force and steric wall- 
repulsion, which lead to migration away from the wall, and a locally varying 
flow-induced orientation, which drives polymer away from the channel center and 
towards the wall. The different dependencies of these effects on the polymer 
bending rigidity and the flow velocity results in a complex dynamical behavior. 
However, a generic effect is the appearance of a maximum in the monomer and 
the center-of-mass distributions, which occurs in the channel center for small flow 
velocities, but moves off-center at higher velocities. 
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1. Introduction 

1.1. Polymers, Vesicles and Cells in Micro- and Nanochannels 

Complex fluids — such as polymer solutions, colloidal dispersions, and suspensions 
of vesicles or cells — exhibit an intriguing flow behavior, particularly in confined 
geometries. The interplay of internal degrees of freedom, e.g, conformational changes 
of polymers or shape changes of vesicles and cells, combined with fluid mediated 
interactions, inhomogeneous flow profiles, and wall interactions lead to novel and 
often unexpected effects [THT5]. Insight into theses aspects provides not only 
an understanding of phenomena for conventional applications such as lubrication, 
adhesion, polymer processing, and blood flow, but also of emergent microfluidic 
devices [19l[20] and the lab-on-a-chip technologies [21] with their micro- or nano-size 
length scale. In this context, the investigation and analysis of single molecules or cells 
is of particular interest, since on the one hand such knowledge is essential to develop 
and optimize micron- and nano-scale devices, while on the other hand the study of 
individual particles often provides a more detailed microscopic knowledge than can be 
extracted from an ensemble. 

The size of red blood cells (RBCs), about S/zm in diameter, is comparable 
to the size of microvessels and microcapillaries in a mammalian body. Thus, the 
deformability and dynamics of an individual cell determines the flow behavior of blood 
in the capillary network. Due to the physiological importance of this process and 
the relatively large size of the objects, detailed investigations of this system started 
already about 50 years ago [22H25] . Modern microfluidic techniques [20|26[27j and new 
mesoscale simulation approaches |28H31| provide new insight into the delicate interplay 
between flow forces, cell deformation, and hydrodynamic interactions [T0 l fr7 l [32 l l33 ] . 
In particular, in microchannels slightly larger than the RBC diameter, it has been 
found that single red blood cells transform from their equilibrium shape of biconcave 
disks into a parachute shape at a critical flow velocity, which depends linearly on 
the bending rigidity and the shear modulus of the cell membrane [10]. At small 
volume fractions of red blood cells, they show a strong clustering tendency due to 
hydrodynamic interactions |17j . At higher volume fractions, several distinct phases are 
found, which range from disordered arrangement of discocytes over single-file motion 
of aligned parachutes to a staggered, zig-zag arrangement of slipper shapes [IT]. In 
channels which are smaller than the RBC diameter, they assume bullet shapes at 
sufficiently high flow velocities [23l|32l[33] . 

Fluid vesicles are distinct from red blood cells since their membrane is fluid, and 
therefore has a vanishing shear modulus. The flow behavior of vesicles has been studied 
intensively in recent years, in particular in shear flows. Much less is known about fluid 
vesicles in capillary flows. In narrow, homogeneous channels, vesicles assume bullet 
shapes [34]. In structured, saw-tooth shaped channels, which are somewhat larger 
than the vesicle size, a complex conformational and dynamical behavior has been 
found [18] . For nearly spherical vesicles, a transition from symmetric shape oscillations 
to orientational oscillations has been predicted with decreasing flow velocity. For 
shapes, which deviate more strongly from a sphere, experiments and simulations show 
new shapes with two symmetric or a single asymmetric tail [18j . 

In narrower channels, with a diameter on the order of 100 nm, it is possible to 
investigate the static and dynamic behavior of single DNA molecules. This system is 
interesting for several reasons. First, microfluidic devices allow for the manipulation. 
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sizing, and sorting of DNA fragments |35[I36| . and thereby a direct visualization of 
genomic information [71[37j- Second, DNA has turned out to be an ideal model 
to study the properties of single polymers, because, when fluoresccntly labelled, 
its conformations, dynamics, and flow properties can be observed directly under a 
microscope [n [38H42] or by correlation spectroscopy [43]. Thus, DNA is also very well 
suited to study the behavior of semiflexible polymers in micro- and nanochannels. 

With a persistence length of Lp ~ llfiui, filamentous actin is a biopolymer for 
which semiflexible behavior becomes important already in wider channels than for 
DNA (with its persistence length Lp ~ 50nm). Thus, experiments with actin are 
essential to elucidate the behavior of semiflexible polymers in microchannels without 
[44H46] and with flow [47]. In equilibrium, the parallel extension of a semiflexible 
polymer shows scaling behavior with power-law dependencies on channel diameter 
and persistence length [JHHHHSI]- 

Structured channels are particular interesting, because they allow to study time- 
dependent flows [T5], threshold forces or flow rates for penetration into narrow 
pores [T61I52H54] . and flow injection into a pore [16]. For fluid vesicles, it has been 
shown that the threshold strength of a driving field (e.g., an electric field) for narrow 
pores increases nearly linearly with membrane bending rigidity and vesicle area, and 
decreases rapidly with the pore radius |52j . For polymers, it has been predicted by 
scaling arguments [53] and confirmed by simulations [TB] that the threshold velocity 
flux for entry into a narrow pore is independent of both the polymer length and the 
pore radius. 

1.2. Polymer Migration in Micro- and Nanochannels 

An interesting aspects of polymer transport in the presence of a zero-slip wall is 
cross-streamline migration. This effect has been observation in planar shear flow 
in the presence of walls as weU as Poiseuille flows [4l|6j|9l[TTl[T2l[T4l[T5l[55] and in 
experiments [SlFflfTS] , as a formation of depletion layer near the wall. The thickness of 
the depletion layer is found to increase with flow, which indicates that migration 
away from the wall increases with the flow rate. This migration phenomenon is 
explained as a results of polymer-wall hydrodynamic interactions. Analytical as well 
as simulation studies, which properly account for such hydrodynamic interactions, 
qualitatively reproduced the experimentally observed wall-induced migration. 

In microchannels, confinement effects and a spatially varying shear rate due 
to the parabolic flow velocity profile imply additional features compared to simple 
shear flows. The steric interaction with the channel wall restricts the conformations 
of a flexible polymer and induce alignment to a semiflexible polymer. Moreover, 
the spatially varying local shear rate changes the conformations and alignment of a 
polymer as a function of its lateral position. In addition to the wall hydrodynamic 
interactions, these two effects also influence the qualitative behavior of cross-streamline 
migration. Simulations of flexible polymers have shown that under strong confinement 
there is a net migration away from the channel center, contrary to the predictions of 
wall induced inward migration jllj . The outward migration is also found in relatively 
weak fiow [15]. In both cases, the outward migration has been attributed to the 
suppression of steric interactions with the wall due to an enhanced alignment by the 
flow. At large fiow strengths, wall hydrodynamic interactions dominate, resulting a 
net inward migration of the polymer center-of-mass. Interestingly, the center-of-mass 
density at the channel center is also found to decrease with flow and at large flow 
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strengths a density maximum is found at a distance away from the center. The reason 
behind the off-center peak of the center-of-mass distribution for flexible polymers is 
the flow-induced conformational change, which leads to stretched polymers close to 
the wall and coiled chains in the central part p[l51f^ im[T^ll41ir5] . This give rise to an 
enhanced outward diffusion and ultimately leads to a concentration dip at the center. 

Although the cross-streamline migration behavior of flexible polymers has been 
extensively studied, there are only a few studies, which focus on the migration behavior 
of semiflexible and rodlike polymers |561l57j. Understanding the migration properties 
of semiflexible and stiff polymers is important for the study of many biopolymers, 
such as actin filaments, microtubules, and intermediate filaments. Since polymer 
conformations are determined by its rigidity, the migration behavior of semiflexible 
and rodlike polymers can differ significantly from that of flexible polymers. 

In this paper, we use computer simulations to study and discuss the migration 
properties of semiflexible polymers with persistence length ranging from half up to 
many times of its contour length. We employ a hybrid simulation scheme including 
Molecular Dynamics (MD) for the polymers and Multiparticlc Collision Dynamics 
(MPC), a particle-based hydrodynamics simulation technique, for the solvent, so 
that wall hydrodynamic interactions and confinement effects are taken into account 
explicitly. The details of the method are given in Sec. [21 Our results, presented 
in Sees. [31 and [H show that polymer migration is indeed qualitatively influenced 
by bending rigidity. Furthermore, in Sec. [5l we investigate the importance of 
hydrodynamic interactions by comparing the results of simulations with MPC and 
Brownian solvents. Finally, the dynamics of polymer migration across the channel is 
studied in Sec. [H 



2. Model and Simulation Method 



2.1. Semiflexible Polymer 



The linear polymer is comprised of Nm point-like monomers of mass M each, which 
are connected by the harmonic potential 

t^^ = y E ~n\~bf, (i) 
1=1 

where ri denotes the position of particle i, b the bond length, and Ks the force 
constant. Excluded- volume interactions are taken into account by the Lennard- Jones 
potential [58] 



Ulj = 



(2) 



(3) 



44(f) -(f) J r<2i/V 
, else 
To account for polymer stiffness, the bending potential is applied [59] 

N^-l 
1=2 

where Ri = — r^-i is the bond vector and is the bending rigidity. In 
the semiflexible limit, the bending rigidity is related to the persistence length by 
Lp = Kb/ksT. 

The dynamics of the polymer is described by Newtons' equations of motion, which 
are integrated using the velocity- Verlet algorithm. 
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Figure 1. Snapshots of polymer conformations close to the channel center (top) 
and the wall (bottom) for the Peclet number Pe = 360 and the persistence lengths 
Lp/Lr ~ 0.5 (a), (d), 2.1 (b), (e), and 30.8 (c), (f). Movies are available as 
Supporting Material I60| . 



2.2. Hydrodynamic Solvent 

In the MPC algorithm, the fluid is described by a set of N point-like particles of mass 
m each, which move in continuous space with velocities determined by a stochastic 
process. Their dynamics evolves in two steps. In the streaming step, the solvent 
particles move ballistically for a time h, which we denote as collision time. In the 
collision step, particles are sorted into the cells of a cubic lattice of lattice constant a 
and their relative velocities, with respect to the center-of-mass velocity of each cell, 
are rotated around a random axis by an angle a. The orientation of the axis is chosen 
independently for every cell and collision step. For every cell, mass, momentum, 
and energy are conserved in this process. The algorithm is described in detail in 
refs. |61H66| . The fluid is confined in a cylindrical channel with periodic boundary 
conditions along the channel axis. No-slip boundary conditions are imposed on the 
channel walls by the bounce-back rule and virtual wall particles, as described in 
ref. |67| . The flow is induced by a gravitational force (mg) acting on every fluid 
particle. 

The interaction of a polymer with the solvent is realized by inclusion of its 
monomers in the MPC collision step [B^. Between two MPC steps, several MD steps 
are performed to update the positions and velocities of the monomers. Extensive 
studies of polymer dynamics confirm the validity of this procedure [58 t [66 l [68 l l69] . 

To maintain a constant temperature, the velocity scaling algorithm is applied 
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as described in Ref. ^70j. Here, an kinetic energy Ei; is chosen from the gamma 
distribution independently for every cell, and the individual particle velocities of 
that cell are multiplied by the factor (2_Efe/(TO ^ w^))^/^, where the sum runs over 
all particles in the considered cell. This assures that the velocity distribution of the 
particles is Maxwellian. Similarly, a Maxwellian velocity distribution is obtained by 
the Monte Carlo procedure described in Ref. [71]. 

Simulations of the pure solvent system yield velocity profiles which agree with 
the solution of Stokes' equation for the considered geometry. 

2.3. Brownian Solvent 

An advantage of the MPC approach is that hydrodynamic interactions can easily 
be switched off, without altering the monomer diffusion significantly |72U73| . In 
this case, denoted as Brownian MPC, each monomer independently performs a 
stochastic collision with a phantom particle with a momentum taken from the Maxwell- 
Boltzmann distribution with variance m (Nc) ksT, where (Nc) is the average number 
of solvent particles per collision cell [73| . 

2.4. Parameters 

For the solvent, we employ the parameters a = 130°, h = O.It, with r = ma^ /ksT 
{ks is Boltzmann's constant and T is temperature), {Nc) = 10, M = m (Nc) ,b = a ^ 
a, the fiuid mass density g = {Nc)m/a^, and ksT/e = 1. The time step in the MD 
simulation is set to Hmd = 5 x IG^-^r. A polymer with = 14 monomers is placed in 
a cylindrical channel of radius R ~ 8a. With the length Lr = {N — l)a = 13a, 
the polymer does not interact with the wall when its center of mass is near the 
channel center. In order to maintain a constant contour length of the polymer, we set 
Kg = 3 X /{ksT/b^). The persistence length Lp is varied, by selecting the values 
Kb/ikBT/b"^) = 7, 28, 50, 100, 200, 400, which correspond to the persistence lengths 
Lp/Lr sa 0.5, 2.1, 3.8, 7.7, 15.4, 30.8. The channel length is 28a. 

The strength of the applied pressure field is characterized by the Peclet number 
Pe = jrji, where 7 ~ ggR/{2r]) is the shear rate at the cylinder wall. The Reynolds 
number Re = gRvm/rj = gR^-j/ {2ri), where Vm is the maximum fluid velocity, depends 
linearly on the shear rate. For the above MPC parameters, the viscosity [BS] is such 
that Re < 1 for all considered 7. Equilibrium simulations for a system with periodic 
boundary conditions yield the end-to-end vector relaxation time th « 3200t. This 
value agrees within approximately 20% with the relaxation time obtained theoretically 
for a semiflexible polymer with the same ratio L,./Lp [74j . The relaxation time of the 
Brownian MPC simulation is t/j w 8200r. 

The averages and probability distributions presented in the following sections are 
calculated in the stationary state for various independent initial conditions. 

3. Radial Distribution Functions 

The imposed flow, with a parabolic flow profile, determines the polymer conformations, 
with a few examples shown in fig.[l] and the polymer distribution functions, such as the 
monomer distribution Pmo and the center-of-mass distribution Pcm- Their dependence 
on flow strength will be characterized in the following sections. 
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Figure 2. Radial monomer distributions for the Peclet numbers Pe = (o), 73 
(□), 147 (o), and 360 (A) of polymers with Lp/Lr 3.8 (a), 7.7 (b), 15.4 (c), and 
30.8 (d). 



3.1. Radial Monomer Distribution 

Radial monomer distributions arc presented in fig. [2] for various flow rates and 
stiffnesses. They are normafized such that 

/ rP„,o{r/a) dr/a^ = 1. (4) 
Jo 

A feature in common with all polymer systems with hydrodynamic interactions is 
the decrease of the monomer concentration adjacent to the wall with increasing flow 
rate. Even without flow, there is a depletion zone for all stiffnesses, which extends 
approximately one radius of gyration into the channel. However, the distributions 
exhibit distinct differences to those obtained previously for flexible and rodlike 
polymers, respectively. For persistence lengths Lp/Lr ^ 8, the density at the center of 
the channel decreases with increasing Peclet number as shown in figs. [2^, b. In contrast, 
for stiffer polymers with Lp/Lj. > 8 the density at the channel center first increases 
with increasing Pe and then decreases at large Peclet numbers. Simultaneously, the 
wall induced migration towards the channel center causes an increase in concentration 
at a finite distance from the center for all stiffnesses. Two effects contribute to the 
formation of the maximum. On the one hand, there is cross-streamline migration 
due to hydrodynamic interactions, and, on the other hand, the flow field causes an 
alignment of a molecule (see section HJ, which increases the local density at the radial 
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Figure 3. Radial monomer distributions for the Peclet numbers Pe = 73 (a) and 
360 (b) and the persistence lengths Lp/Lr ^ 0.5 (+), 2.1 (•), 3.8 (□), 7.7 (o), 
15.4 (A), and 30.8 (*). 



positions of strong flow-induced alignment. 

The distribution functions are rather similar for the various stiffnesses at Pe = 
within the accuracy of the simulation. Interestingly, stiffer polymers exhibit a more 
pronounced migration away from the wall for Pe < 100 as is shown in fig. [3^. 
Simultaneously, Pmo decreases in the channel center for small stiffnesses, leading to an 
off-center maximum, and increases for larger ones, where the maximum is at the center 
for some of the larger stiffnesses. With increasing flow rate, off-center density maxima 
are obtained for all stiffnesses. For the flow rate Pe = 360, wc observe a decrease in 
the depletion zone adjacent to the wall with increasing stiffness for Lp/Lr < 8 and a 
reversion of the trend at larger stiffnesses as depicted in flg. [31d. 

The width of the distribution, 

(r^) = / r^'Prnoir/a) dr/a\ (5) 
Jo 

emphasizes the differences in migration behavior. As displayed in fig. [6^. stiffer 
polymers exhibit a smaller width at moderate Peclet numbers, indicating an enhanced 
net inward migration with increasing bending rigidity, but (r^) seems to saturate at 
large flow rates. In contrast, the widths of more flexible polymers show only minor 
changes at small Pe and decrease more rapidly at larger Peclet numbers. Their widths 
can be smaller than those of stiffer polymers. 

3.2. Radial Center-of-Mass Distribution 

Radial polymer center-of-mass distributions arc displayed in flg. |4l Qualitatively, Pcm 
exhibits similar features as the monomer distribution Pmo, reflecting the same physical 
mechanisms. For Lp/L,. < 8, the concentration at the channel center decreases with 
increasing Pe. The density profiles for 50 < Pe < 150 are very similar and indicate 
a weak dependence of the center-of-mass properties on the flow rate only. For high 
Pe values, the densities at the channel center decrease signiflcantly for all stiffnesses, 
and a clear maximum appears at rem > 0. The density decrease with flow near the 
wall is mainly due to wall induced cross-streamline migration of polymers and partly 
due to steric interactions of the very stiff polymers with the wall, while the density 
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Figure 4. Radial ccntcr-of-mass distributions for the Pcclct numbers Pe = (o), 
73 (□), 147 (o), and 360 (A) of polymers with Lp/Lr « 3.8 (a), 7.7 (b), 15.4 (c), 
and 30.8 (d). 





Figure 5. Radial center-of-mass distributions for the Peclet numbers Pe = 73 
(a) and 360 (b) and persistence lengths Lp/Lr ^ 0.5 (+), 2.1 (•), 3.8 (□), 7.7 (o), 
15.4 (A), and 30.8 {*). 
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Figure 6. Widths (r^) of the radial monomer distributions (a) and centcr-of- 
mass distributions (b) for the persistence lengths Lp/Lr ~ 0.5 (•), 2.1 (□), 3.8 
(o), 7.7 (A), 15.4 W, and 30.8 {+). 



decrease at the channel center indicates a migration away from the center. The steric 
contribution to migration is evident from the Brownian MPC simulations presented 
in sec. [5l The two competing migration mechanisms result in the formation of a 
maximum at a distance away from the center at large Pe. The maximum is less 
pronounced for the center-of-mass distribution than for the monomer distribution. 

Figure [5] shows center-of-mass distributions of polymers with different bending 
rigidities for Pe « 73 and Pe « 360. As for even smaller Peclet numbers, the 
distributions exhibit a maximum essentially in the channel center for Pe « 73. An off- 
center maximum only appears for larger flow rates, which is accompanied by a density 
decrease in the channel center. Moreover, a pronounced depletion layer is visible at 
the wall for large stiffnesses. 

Flexible polymers exhibit a pronounced density increase adjacent to the wall with 
increasing flow rate and before migration sets in, whereas very stiff polymers exhibit 
an increased depletion layer at the wall with increasing flow rate, which we attribute 
to steric polymer-wall interactions. This is reflected in the flow rate dependence of 
the width of the center-of-mass distribution displayed in fig. [6j). For Lp/Lr ^ 4, 
the width increases with increasing Peclet numbers for small Pe. This increase is 
explained by the suppression of polymer-wall steric interactions due to flow alignment 
of polymers (cf. section HJ . For flexible polymers such an alignment is enhanced 
by flow-induced polymer stretching along the channel axis. However, at larger Pe, 
the wall induced hydrodynamic lift force dominates the dynamics, which results in 
a decrease of width. This dependence is qualitatively different for polymers with 
Lp/Lr ^ 5, where wall induced migration dominates even at low Pe, implying a 
decrease in width. However, at large Peclet numbers, we observe a saturation or even 
increase in width by outward migration due to hydrodynamic interactions, which 
counter balances the radially inward migration, as discussed in Sec. 6. 

4. Alignment and Conformations 

4.1. Orientational Order Parameter 

In microchannel flows with no-slip boundary conditions the local shear rate changes 
linearly with radial position. Hence, the flow-induced force experienced by a polymer 
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depends on its radial position. Its conformations and alignment, in response to the 
flow forces, is expected to vary with its bending rigidity and its radial center-of-mass 
position. Studies of flexible and semiflexible polymers [T51I75] reveal large orientational 
changes by the imposed flow, which is important for cross-streamline migration |55j . 
To study the influence of flow on the polymer orientation, we consider the orientational 
order parameter 

Sir,^,) = ^{3 cos-" 0-1), (6) 

where 6 is the angle between the polymer end-to-end vector and the flow direction, 
as a function of its center-of-mass radial position. S{rcm) is plotted for various 
Peclet numbers and stiffnesses in fig. [71 In the absence of flow and for all Kf,, the 
polymer orientation is isotropic at the channel center. For distances close to the 
wall, confinement causes a preferred orientation along the channel axis with 5' ~ 1. 
With increasing flow rate, the order parameter increases for all stiffness and reaches 
a plateau at large Pe for certain radii. The plateau value and its extension depends 
on stiffness as shown in flg. [51 

The order parameter exhibits a non-monotonic dependence on flow strength close 
to the channel center and depends on polymer stiffness as shown in the inset of flg. [51 
An increase in stiffness leads to a decrease of the order parameter for Lp/Lr < 15. For 
larger stiffnesses, S{0) is larger for larger Lp/Lr at aU radial distances. At Pe = 360, S 
decreases with increases stiffness and becomes even negative for Lp/Lr — 2, compare 
the inset of fig. [51 As the polymer stiffness increases, S increases close to the channel 
center. The negative values of S at the channel center are attributed to the formation 
of U-shaped conformations, compare flg. [H Such conformations are not possible for 
flexible polymers |15| and polymers with large bending rigidity, hence S assumes higher 
values at the center as Kb increases. 

4-2. Bending Energy 

U-shaped conformations can be quantitatively analyzed by calculating the average 
bending energy Ub- Figure [9l displays C/f, as function of the radial polymer center- 
of-mass position for various persistence lengths and for various Peclet numbers. For 
Pe = the bending energy is nearly uniform across the channel cross section and 
decreases near the wall, due to wall induced alignment. The magnitude of Ub is close 
to the thermal average Ub = {Lp/a — l)kBT of the nearly harmonic bending potential 
for Lp/Lr > 1, as expected. An increase in Pe results in an increase in its absolute 
value, which is a result of the conformational changes of the polymer in response to 
the force exerted by the flow, compare flg. [1] At a given Pe and radial position, the 
flow exerts the same force on a polymer, but its conformations depend on its bending 
rigidity. Polymers with large bending rigidity hardly undergo conformational changes, 
while polymers with rather small Kb go through signiflcant conformational changes, 
but their bending energies are similar. This is exemplifled in fig. I10[ where Ub/iksT) 
is shown for various bending rigidities and Pe — 360. 

In figs. [HI and [TUl the radial energy profiles are qualitatively similar for all Kb 
(with Lp/Lr > 1) for a given Pe. For 0.5 < Lp/Lr £ 5, the polymer adopts U- 
shaped conformations at the center and hence Ub is large. Away from the center, 
Ub decreases as U-shapes disappear. With increasing radial center-of-mass position, 
the local shear rate increases and a polymer with small persistence length assumes 
transient bent conformations. The average over individual configurations provides a 
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Figure 7. Orientation order parameter for tlie Peclet numbers Pe = (o), 73 
(□), 147 (o), and 360 (A) of polymers witli L^/Lr 3.8 (a), 7.7 (b), 15.4 (c), and 
30.8 (d). 
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Figure 8. Orientational order parameters S for Peclet numbers Pe = 73 (a) and 
360 (b) and persistence lengths Lp/L^ ^ 0.5 (+), 2.1 (•), 3.8 (□), 7.7 (o), 15.4 
(a), and 30.8 (*). Inset: S close to the channel center for Pe = 360. The line is 
a guide to the eye. 
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Figure 9. Bending energies Ut, for Peclet numbers Pe = (o), 73 (□), 147 (o), 
and 360 (A) of polymers with persistence lengths Lp/Lr ~ 3.8 (a) and 30.8 (b). 
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Figure 10. Bending energies for Peclet number Pe = 360 and persistence lengths 
Lp/Lr ^ 0.5 (•), 2.1 (□), 3.8 (o), 7.7 (a), 15.4 (*), and 30.8 (+). 



high value for Ub for such rem- The extent of bending decreases as Kb is increased, since 
for stiffer polymers (1 — cosi?) ^ Ub/Kb, where "(9 is the angle between successive bond 
vectors. More aligned bonds and hence smaller angle i9 imply that Ub/Kb decreases 
with increasing Kb, in agreement with the results of fig. IHl 



5. Structural Properties w^ithout Hydrodynamic Interactions 

As is generally accepted by now, hydrodynamic interactions determine the behavior 
of polymers in microcapillary flows. In order to separate effects due to intramolecular 
and surface hydrodynamic interactions from those caused by the flow and steric 
interactions, we performed Brownian MPC simulations [73| . Brownian dynamics 
simulation have also been performed in Ref. |57) , with a polymer composed of rodlike 
segments. However, intramolecular hydrodynamic interactions are taken into account 
by employing anisotropic drag coefficients parallel and perpendicular to the rodlike 
segments. In contrast, our polymer model is composed of point-like monomers with 
isotropic friction. Hence, we consider a free draining chain, whereas in Ref. |57j 
intramolecular hydrodynamic interactions are taken into account implicitly. 
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Figure 11. Radial monomer distributions for systems without hydrodynamic 
interactions. The persistence lengths are Lp/Lr = 3.8 (a) and 30.8 (b), and the 
Peclet numbers Pe = (o), 147 (□), 360 (a), and 920 (*). 
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Figure 12. Radial ccntcr-of mass distributions for systems without 
hydrodynamic interactions. The persistence lengths are Lp/Lr = 3.8 (a) and 
30.8 (b) and the Peclet numbers Pe = (o), 147 (□), 360 (a), and 920 (*). 



5.1. Radial Concentration Distributions 

The radial monomer probability distributions for the persistence lengths Lp/Lr ~ 2.1 
and 30.8 are presented in fig. [11] for various flow rates. For Pe < 200, the distribution 
functions are rather similar. Adjacent to the wall, the distributions for small stiffnesses 
are even similar for all flow rates. At large flow rates and persistence lengths 
Lp/Lr < 3.5, we find a monomer density increase at r^m/ o, ~ 4, which we attribute to 
the induced alignment of the polymers by the fiow. However, for larger stiffnesses, we 
find an inward migration adjacent to the wall due to increased steric polymer-surface 
interactions (cf. flg. [Tib). Here, the polymer tumbling motion combined with the 
smaller number of conformational degrees of freedom at large stiffnesses leads to a 
strong repulsion from the wall. At the same time, the density in the channel center 
increases with a maximum at r w 0. 

The radial center-of-mass distribution exhibits an outward migration even without 
hydrodynamic interactions for Lp/Lr < 3.5 due to polymer alignment (cf. flg. [T2|). In 
contrast, Pcm becomes narrower with increasing fiow rate for Lp/Lr ~ 30.8, similar 
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to the monomer distribution. Hence, a qualitative different behavior is obtained 
depending on the abiUty of the polymer to adjust to the conformational restrictions 
by the walls. 

The Brownian simulations of ref. |57| display in strong outward migration of the 
polymers, which leads to a density minimum at the channel center and a large density 
increase near the walls. Two effects contribute to this behavior. On the one hand 
alignment of the polymers and on the other hand anisotropic diffusion by the aligned 
segments. We do not observe such an increase in density adjacent to the wall, because 
no hydrodynamic effects are included in our simulations. 

5.2. Alignment 

The orientational order parameters display qualitatively similar dependencies on flow 
rate as those for systems with hydrodynamic interactions (cf. figs. [7] and [5]). However, 
alignment is far less pronounced and much larger Peclet numbers are required to 
achieve significant flow alignment. For Pe = 360, we find S w 0.2 at Lp/L,. « 3.8, 
compared to S = 0.6 in the presence of hydrodynamic interactions. At the same time, 
we also observe the appearance of U-shaped structures in the channel center for larger 
Peclet numbers. The comparison shows that certain qualitative features of polymer 
alignment are determined by the flow profile, rather than hydrodynamic interactions. 
However, hydrodynamic interactions clearly influence the polymer orientation in a 
quantitative manner and lead to a more pronounced alignment. 

6. Cross-Streamline Dynamics 

The most striking effect of flow on a semiflexible polymer is the strong dependence of its 
orientation on Pe and the radial distance. This aspect provides the key to understand 
the appearance of the off-center maximum in the center-of-mass distribution function. 

As is well known, a rod in solution exhibits a larger diffusion coefficient parallel to 
its axis than perpendicular to it. In the absence of flow, a semiflexible polymer behaves 
very similar to a stiff rod, we therefore expect that the observed polymer orientational 
differences across the channel will lead to differences in the lateral diffusion behavior. 

In order to characterize this diffusive behavior, we calculate the radial 
mean displacement Ar(rcm) = ([rcm(to + At) - rcr„(to)] ecm(io)), where Scm, = 
fcm,/\fc'm\) a-nd the center-of-mass mean square displacement Ar^(rt,m) = 
{(Tcmito + At)ecm(to) - (^■cm(io)ecm(io)))^) within a certain time At for a polymer 
with the center-of-mass position rem at to- We choose the time At such that the 
center-of-mass displacement is smaller than a and the mean square displacement is 
linear in time. For an homogeneous and isotropic system, Ar and Ar^ are independent 
of the initial position. 

The mean values are the first and second moment of the conditional probability 
distribution function to find the polymer center-of-mass at position rcm{to + At) at 
the time tg + At, when it has been at position rcm{to) at the time to. Initially, 
the distribution is very narrow, ideally a 5 function, and broadens in time. In the 
limit At — >■ cx) the conditional probability distribution turns into the equilibrium 
distribution function as displayed in figs. |4] and [H 

Figures [T3l and [T4l displav Ar and Ar^ for various Peclet numbers and stiffnesses. 
The mean displacements evidently depend on flow rate and stiffness. At Pe = 0, 
there is no drift in the channel center. Only for rcm/o, > 1.5, we observe a drift, which 
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Figure 13. Radial centcr-of-mass displacements Ar(r'cm) and mean square 
displacements Ar^{rcm) (insets) for Peclet numbers Pe = (o), 73 (□), 147 
(o), and 360 (a) of polymers with persistence lengths Lp/Lr ~ 3.8 (a) and 30.8 
(b). 




Figure 14. Radial center-of-mass displacements Ar{rcm) and mean square 
displacements Ar^(rcm.) (insets) for persistence lengths Lp/Lr ~ 0.5 (-I-), 2.1 
(•), 3.8 (□), 7.7 (o), 15.4 (A), and 30.8 {*), and Pe = 73 (a) and Pe = 360 (b). 



we attribute to steric polymer wall interactions, and which increase adjacent to the 
wall. For larger flow rates, the drift for rem ^ 4 is close to zero due to alignment of a 
polymer and lack of steric wall interactions (cf. fig. [14] a). At large Peclet numbers, 
Ar is positive for rem/ a < 2.5 (cf. fig. [Mb) and all persistence lengths Lp/Lr > 2, 
which indicates an outward migration. The wall interactions lead to a larger drift 
of polymers towards the channel center adjacent to the wall, which increases with 
increasing stiffness. 

The mean square displacement decreases with increasing radial distance for 
fern < 4a at all stiffnesses and flow rates. Since Ar is rather small for r^m < 4a, 
Ar^ is determined by diffusion, i.e., the radial diffusion is larger in the channel 
center. This leads to a minimum in the center-of-mass distribution at large flow 
rates and stiffnesses. At larger radial distances and large stiffnesses, the mean 
square displacement is strongly affected by wall interactions. There is no simple 
explanation for the non-monotonic behavior of the mean square displacement. Various 
wall interactions contribute to the system behavior: the wall lift force, steric wall 
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Figure 15. Radial center-of-mass mean square displacements without 
hydrodynamic interactions for the persistence lengths Lp/Lr ~ 30.8 and the 
Peclet numbers Pe = 360 (A), 920 (*), and 1830 (x). 



interactions, and an increase of the tumbling frequency with increasing flow rate 



Without hydrodynamic interactions, the drift at Pe = is comparable to that 
of the system with hydrodynamic interactions. It changes with increasing Peclet 
number due to flow alignment, but only slightly. The mean square displacement 
is independent of the radial position for distances rem < 6a for all Peclet numbers 
Pe < 1000 and stiffnesses as shown in fig. [121 Only close to the wall, there is a decrease 
due to steric wall interactions. This supports our conclusion that intramolecular 
hydrodynamic interactions are responsible for outward migration in systems with 
hydrodynamic interactions. At much larger Peclet numbers (Pe ~ 2000), the mean 
square displacement is constant for r < 4a. increases then with the radial distance 
and decreases again beyond r = 6a, similar to the behavior shown in the inset of 
fig. [Hk for Lp/Lr ~ 15.4. This behavior is explained by the strong steric polymer 
wall interaction and fast tumbling motion at high fiow rates. Note that the tumbling 
frequency increases with the shear rate |741I75| . Hence, systems with and without 
hydrodynamic interactions exhibit similar features, but larger Peclet numbers are 
required for systems without hydrodynamic interactions and the behavior of systems 
with hydrodynamic interactions is richer and additional features are displayed. 

7. Summary and Conclusions 

We have analyzed the fiow behavior of semiflexible polymers confined in a cylindrical 
channel and systematically studied the influence of bending rigidity on structural 
properties and cross-streamline migration. The study includes polymers with 
moderate bending rigidity, comparable to actin fllaments, to very stiff rodlike 
polymers. 

Poiseuille flow induces a strong radial-dependent polymer alignment parallel to 
the channel associated with a radially outward migration of its center-of-masses at 
small Peclet numbers and, at large Peclet numbers, a wall induces inward migration. 

Stiffness changes the polymer behavior qualitatively and quantitatively. Taking 
a system without hydrodynamic interactions as a reference, there is no outward 
migration and there arc no wall hydrodynamic forces. Correspondingly, the center- 
of-mass density decreases in the channel center with increasing flow rate and up to 
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a certain stiffness {Lp/L,. < 20) due to polymer alignment. For large stiffnesses, 
however, steric polymer-wall interactions cause inward migration of the molecule and 
an increase in density at the channel center with increasing flow rate. We attribute 
this effect to the increase of the tumbling frequency at large flow rates and hence more 
polymer-wall contacts. Hence, even for such a system, a qualitative different behavior 
is found for large and small stiffnesses. 

This picture is enriched by hydrodynamic interactions, which enhance outward 
migration. At the same time wall-lift forces lead to inward migration and an off- 
center density maximum is formed at large flow rates. Such a maximum has also 
been found in experiments using actin filaments |47j . The center-of-mass distribution 
strongly depends on stiffness and the presence of hydrodynamic interactions, as is 
reflected in its width. Flexible polymers exhibit a more complex migration behavior 
than stiffer polymers due to conformational changes, i.e, stretching and alignment. 
For large stiffnesses, wall steric interactions lead to a strong inward migration as 
for polymers without hydrodynamic interactions. The semiflexible polymers exhibit 
a non-universal behavior when changing stiffness and flow rate. At low flow rates, 
stiffer polymers display a stronger inward migration, whereas at high flow rates a 
reversed trend is observed initially and ultimately a strong inward migration at large 
stiffnesses. It is the interplay between confinement, hydrodynamic interactions, and 
stiffness which leads to the emergence of a complex structural and dynamical behavior 
of the polymer. This is reflected in the stiffness dependence of the mean radial and 
mean square radial displacements. 

Although the understanding of the flow behavior of soft mesoscale particles in 
nano- and microchannels has made considerable progress in recent years, this field 
can be expected to attract far more attention in the future. For example, flow in 
structured channels has hardly been addressed so far. Also, interactions between 
polymers at flnite concentration, and between polymers and other soft particles, needs 
to be investigated. 
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